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Abstract 

We report preparation, x-ray diffraction, magnetic susceptibility x{T) ^-iid heat capacity Cp{T) 
measurements on the undoped samples as also samples with Zn-doped (5 = 0) at Cu site, Ni- 
doped {S = 1) at Cu site, and Ca-doped (holes) at Bi site in the coupled two-leg spin ladder 
system BiCu2P06. While, Zn shows complete solid solubility, Ni could be doped to about 20% and 
Ca to about 15%. Magnetisation and heat capacity data in the undoped compound point towards 
the existence of frustration effects. In all the samples, the xiT) at low temperature increases 
with doping content. The Zn- induced susceptibility is smaller than that due to effective 5 = 1/2 
moments possibly due to frustrating next-nearest-neighbour interactions along the leg. For Zn 
content x > 0.01, xiT) deviates from the Curie-law at low temperatures. The magnetic specific 
heat data Cm{T) for the Zn-doped samples show weak anomalies at low temperature in agreement 
with x{T) behavior. The anomalies are suggestive of spin freezing at low-T. In contrast, prominent 
effects are observed in x{T) and Cm{T) on Ni-doped samples. The zero-field-cooled (ZFC) and 
field-cooled (FC) x(^) data are different from each other at low temperature unlike that for Zn- 
doped samples, clearly indicating a transition to a spin-glass like phase. No anomalies were found 
in Ca- or Pb-doped samples. 

PACS numbers: 75.40. Cx, 75.10.Pq, 75.10.Hr 
Keywords: Spin ladder, BiCu2P06, spin-glass 
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I. INTRODUCTION 

Currently, there is great interest in the magnetic properties of pure and doped low- 
dimensional Heisenberg antiferromagnetic (HAF) and frustrated spin systems. Exotic 
ground states such as valence bond solid (VBS) state in Haldane chains, spin-dimers, and 
spin-Peierls systems^ or the resonating valence bond (RVB) state in chain, ladder^ or planar 
systems^ have been investigated. A lot of experimental and theoretical work has been done 
on the effect of non-magnetic Zn^"*" {S — 0) and magnetic Ni^+ {S — 1) impurities on Cu- 
based systems. When Cu^"*" ions in hole-doped, metaUic Cu02 planes in YBa2Cu306+5 are 
substituted with Zn^+, nuclear magnetic resonance (NMR) experiments^ revealed that local 
magnetic moments are induced on neighboring Cu sites. Further NMR studies^ investigated 
the effect of impurities on the magnetic correlations and showed that each vacancy causes 
a strong antiferromagnetic polarization around its immediate vicinity.^ This phenomenon 
was also observed in Ni-doped square lattices.^ Similar behavior was also found in Haldane 
chains^ and ladders. 

In the case of 5" = |, even-leg spin ladders (which have a spin gap), a transition from 
the spin-liquid ground state to an ordered state across a quantum critical point (QCP) can 
occur as a function of some tuning parameter such as the interladder coupling.^ Additionally, 
non-magnetic impurities have been found to induce disordered magnetism or magnetic long- 
range-order (LRO) in a prototypical two-leg ladder such as SrCu203.^^'^^ In SrCu203, the 
exchange interaction along the leg (Ji/Zcb) was found to be 1900 K. The ratio of the rung 
exchange to leg exchange J2/J1 and the spin-gap (A/A;^) were found to be 0.5 and 450 K, 
respectively, with only a weak inter-ladder interaction.^^ In this spin gapped system, even 
a small amount of disorder by introduction of non-magnetic Zn^"*" {S = 0) impurities at 
the Cu^+ site induces long-range AF ordering at the Neel temperature Tn ranging from 3 
to 8 K for 0.01 < a; < 0.08.^^ The specific heat measurements above for x = 0.02 and 
X — 0.04 also confirmed the gapless antiferromagnetic state from its linear variation with 
temperature, which should be observed in gapless systems such as uniform spin chains. 
Similar AF ordering was also found in Ni doped SrCu203 from NMR experiments. ^° On the 
other hand, in the systems which are beyond the QCP, i.e., on the LRO side of the phase 
diagram such as the 2D square planar La2Cu04, doping with Zn leads to a reduction of T^.^^ 
In contrast, the physics of defect induced moments in geometrically frustrated spin systems. 
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such as kagome, '^^ hyper-kagome systems^^, etc. is quite different from non-frustrated spin 
systems. These interesting results inspired us to look for new quantum spin systems and 
investigate the impurity effects in them. 

Recently, we have reported the magnetic susceptibility and the specific heat of a new 
spin-ladder system BiCu2P06, in which the exchange interaction along the leg Ji/ks was 
found to be 80 K and the ratio of the rung to leg exchange coupling J2/J1 to be nearly 
1.^° Prom these exchange couplings and magnetic specific heat, we estimated a spin-gap 
(A/Zcs) of 34 K, which is in good agreement with previously reported value of 33 K.^^ The 
first-principles band structure and magnetic susceptibility analysis^° suggested a strong inter- 
ladder coupling J3/J1 ~ 0.75 in the 6c— plane and a next-nearest-neighbor (nnn) exchange 
interaction along the leg J4/ Ji ~ 0.34. While we are not aware of any calculations of the QCP 
for 2D coupled two-leg ladders with a frustrating nnn interaction along the leg, J3/ Ji ~ 0.75 
would seem to be very large and presumably close to a QCP. Perturbing this ground state by 
doping with non-magnetic or magnetic impurities can reveal important information about 
the nature of the system. 

Furthermore in an earlier paper, we have used local probe measurements (nuclear mag- 
netic resonance NMR and muon spin resonance /iSR) to probe impurity (Zn and Ni) cfi^ccts 
in the BiCu2P06 spin-ladder. We have shown the appearance of impurity-induced moments 
leading eventually to spin-freezing at low temperatures. We now report bulk (magnetic sus- 
ceptibility and heat capacity) measurements on pure, and doped (homovalent Zn, Ni at Cu 
site and heterovalent Pb, Ca at Bi site) HAF spin-ladder BiCu2P06 i.e. Bi(Cui_j;Zna;)2P06, 
Bi(Cui_.yNij,)2P06 and Bii_2Ca2Cu2P06. While such measurements vahdate the sample 
quality and surely constitute the first steps to investigate the gross properties of a system 
such as the occurence of a spin-gap, phase transitions, etc., they also enable us to reach 
conclusions about details such as the applicability (or lack thereof) of theoretical mod- 
els (isolated spin-ladder, Majumdar-Ghosh chain, dimer model, etc.) and the value of the 
impurity- induced magnetic moment. These are useful inputs to theorists to refine and 
improve their theoretical models. Our magnetic susceptibility measurements on Zn-doped 
BiCu2P06 evidence low-tempcraturc upturns due to induced magnetic moments. The data 
also suggest that the Zn-induced moment is smaller than that expected for isolated lad- 
ders and is likely related to frustrations present in our system. In the case of Ni doping, a 
transition to a spin-glass-like disordered state is clearly evident in the susceptibility data. 
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Anomalies were also observed in the heat capacity data. Similar to Zn-doped samples, in 
Ni-doped samples the induced moment (combination of the Ni magnetic moment and the 
neighbouring staggered moments) is found to be smaller than that for an isolated -S" = 1 
entity. No anomalies were found in the susceptibility or heat capacity data of Ca or Pb 
doped samples. 

II. EXPERIMENTAL DETAILS 

Polycrystalline samples of BiCu2(i_j;)Zn2j;P06, {0 < x < 1), BiCu2(i_y)Ni2yP06, (0 < y < 
0.2), Bii_^Pb^Cu2P06, (0 < ^ < 0.07) and Bii_^Ca^Cu2P06, {0 < z < 0.15) were prepared 
by conventional solid state reaction methods, using Bi203, CuO, ZnO, NiO, CaCOs and 
(NH4)2HP04. Samples were ground, pelletized, and fired at 770°C for four days with three 
intermediate grindings. X-ray diffraction measurements were performed using a PANalytical 
X'pert PRO powder diffractometer. All the above samples were found to be single-phase 
and structural parameters were refined using FuUprof software.^^ Magnetization M was 
measured as a function applied magnetic field H and temperature T {2 < T < 300 K) 
using a SQUID magnetometer and a vibrating sample magnetometer (VSM) of a Physical 
Property Measurement System (PPMS), both from Quantum Design, Inc. The specific 
heat Cp measurements were performed as a function of T in the range 1.8 to 70 K by the 
relaxation method using the PPMS. 

III. RESULTS AND ANALYSIS 

A. Crystal structure and lattice parameters 

In the crystal structure of BiCu2P06,^^ two edge-shared CuO^ distorted square pyramids, 
which are built from two different Cu (Cul and Cu2) and O (01 and 02) atoms, respectively, 
form zigzag two-leg ladders along the crystallographic 6— direction. The bond angle and bond 
length along the leg (Cul-01-Cu2) of the ladder are about 112° and 3.2 A, respectively and 
along the rung (Cul-02-Cu2) the corresponding values are 92° and 2.8 A, respectively. The 
Bi atoms are positioned between the ladders in the 6c— plane as shown in Fig. 1(a). These 
planes are separated by PO4 tetrahedron units. Various exchange couplings {Ji, J2, J3, and 
J4) between Cu atoms in its fee-plane have been indicated in Fig. 1(b). 
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All the above mentioned samples were found to be single phase from x-ray diffraction 
(xrd) measurements and Riet veldt refinement analysis. The Riet veldt refinement pattern 
for the undoped sample is shown in Fig. 2(a). The xrd patterns for Zn doped samples were 
also indexed based on the orthorhombic space group Pnnia. The xrd peaks shift towards 
smaller angles 26 as the Zn content changes from a; = 0.0 to x = 1.0 which is a clear 
indication of a change in the lattice parameters. The yielded lattice parameters a, b, and 
c from the refinement are plotted as a function of x in Fig. 2(c). The lattice parameters 
of BiCu2P06 {x — 0) and BiZn2P06 {x — 1) are in reasonable agreement with previously 
reported data^^ and the values for intermediate concentrations also smoothly interpolate 
between the end points. The changes in lattice parameters on Ni, Pb, and Ca doping (not 
shown here) are somewhat less than that on Zn doping . 

IV. MAGNETIC SUSCEPTIBILITY AND HEAT CAPACITY 
A. Undoped BiCu2P06 

The temperature dependence of bulk magnetic susceptibility x(T') = M/H for BiCu2P06 
is shown in Fig. 3. The xiT) shows a broad maximum at about 56 K, indicative of a low- 
dimensional magnetic system, thereafter decreasing steeply with decreasing temperature. 
This is indicative of a gap in the spin excitation spectrum in agreement with the 2-leg 
ladder geometry suggested by the structure. A small upturn appears below 7 K, likely due 
to extrinsic paramagnetic impurities and /or natural chain breaks in our polycrystalline 
samples. We first fit the low-temperature part of x(T) below 15 K to the equation^^ 



This yields the parameters Xo = (4.4 ±0.2) x 10"^ cm^ /mole Cu, C = (3.6 ±0.2) x 10"^ 
cm=^ K /mole Cu, ^ ~ K, ^ = ^f^'^ = (2.6 ± 0.1) x 10"^ and the spin-gap A/ke = 



(45±3) K. Here is the Boltzmann constant, is the Avogadro number, g is the Lande'g'- 
f actor (fixed to 2.1, a typical value in cuprates), and 7 is the curvature of magnon dispersion 
at the band minimum. The temperature independent susceptibility Xo is the sum of the core 
diamagnetic susceptibility Xcore, the Van Vleck susceptibility Xw ^-nd a residual susceptibility 




(1) 
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Xres (this will be explained shortly). The Xcore is calculated^^ to be —0.6 x 10~^ cm^/mol 
Cu. Our recent NMR results^^ have shown that the ^^P NMR shift (and hence the spin 
susceptibility) falls exponentially at low-T due to the gap, there remains a non-zero NMR 
shift as T ^ 0. This amounts to a residual spin susceptibility Xres = 1-4 x 10~"^ cm"^/mol 
Cu. Using this value of Xres we get Xw = 3.6 x 10~^ cm^/mol Cu which is somewhat higher 
than Xvv of other cuprates. The Curie constant corresponds to about 0.1 % of isolated spin- 
I impurities. In a latter section, we analyse the impurity-induced susceptibility in greater 
detail. The ratio 7/ A was also found from the fit. However, due to the small amount of 
data at low-T, the result was found to be dependent on the T-range used for the fit and is 
therefore considered unreliable. 

In order to gain further insight into the magnetic interactions present in this system, 
we now attempt to fit the susceptibility data in the full temperature range for undoped 
BiCu2P06. Note that while the measured susceptibility has a low-T Curie term (albeit 
small) the spin susceptibility obtained via the ^^P NMR shift measurements {xnmr{T), 
detailed in a separate paper^^) is intrinsic and has no Curie term. We have, therefore, used 
Xnmr (T) to carry out the analysis following Equation (2) below.^^ This kind of equation 
has been used in the past to determine the strength of inter-dimer, inter-chain, and inter- 
ladder interactions in the corresponding systems (non-frustrated) like SrCu203,^'' CaV205,^^, 



Here, xJaddcr ('^i,-'2, 7") = ^^pTIJ" gladder (T) (with Xiadder(T) being the isolated ladder 
susceptibility^*) while A (= 2{zeff — Zq)) accounts for the intersystem exchange coupling. 
The fitting parameters depend somewhat on the T-range used for the fit. The obtained 
values for the leg and rung coupling are about Ji/ks ~ 80 K & J2/kB — 78 K. Also, A is 
found to be in the range 8 — 10. Given that 2^0, the magnetic coordination number, in the case 
of the isolated ladder is 3, we get the effective magnetic coordination number z^jf ~ 7 — 8. 
Using this and our experimental value of t^"*"^ = Xnmr — 1-25 x 10~^ cm^/mole Cu, we 

get iyg2^i " = 0.42 — 0.48 which lies a little above the universal curve (for nonfrustratcd 

bipartite AF spin lattices) given by Johnston (see Fig. 4 of Ref.^^). Further, Johnston^^ 

-.max J max ^ 

had made the observation that - — 2 -2 hes above the universal curve for highly frus- 



SraCrsOg,^^ etc. 
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trated systems like the triangular system and other cluster systems where the frustrating 
interaction is equal to the other interaction. One might, therefore, naively conclude that 
BiCu2P06 might be only moderately frustrated. 

The plot of XspinT (= (x ~ Xcore — Xvv)T) vs T for undopcd BiCu2P06 is shown in Fig. 
3(b). Even at 300 K, which is nearly four times the value of J^^^/ks, the value of XspinT 
is about 0.185 and smaller than 0.375, which is the value of XspinT in the paramagnetic 
limit. This suggests that the system is not in the paramagnetic hmit even at 300 K. As a 
further illustration, we show the spin magnetization Mspin{H) as a function of applied field 
at T = 350 K with the comparison of Brillouin function for 5" = | and g — 2.1 in the inset 
of Fig. 3(b). Clearly, the data arc well below the Brillouin function curve. These results 
suggest that AF and frustrating interactions are present even at T = 350 K. 

Theoretically, a coupled two-leg spin ladder system (Ji = J2) has a quantum critical 
point (QCP) at about J3/J1 = 0.3, which separates the spin-gap region from the LRO 
antiferromagnetic region.^ Experimentally, however, we have proven that this system has a 
spin gap of about 35 — 40 K. Recent inelastic neutron scattering results on BiCu2P06 yield 
a spin-gap of about 1.5 meV (equivalent to 20 K).^'^ The difference with our results could 
be due to a wave- vector dependence of the spin-gap affecting the neutron scattering results 
on powder samples. Prom theory it is known that in two-leg ladders the existing spin-gap 
increases due to frustrating nnn interactions.^^ Also, from theoretical DMRG calculations, 
it is predicted that the three-leg spin ladder (which is gapless) will develop a spin-gap 
in case of frustrating nnn interactions. The continued presence of the spin-gap state in 
our system inspite of the relatively large inter-ladder couplings (predicted from electronic 
structure calculations^'^) points to the presence of frustrating interactions. It is hkely that 
the next-nearest-neighbour {nnn) coupling J4 is frustrating and plays a role in sustaining 
the spin-gap state. 

Specific heat Cp{T) of BiCu2P06 is shown in Fig. 4. We have analyzed Cp{T) of 
BiZn2P06, the non-magnetic analog of BiCu2P06, using the Debye model. The data 
couldn't be fit using a single Debye temperature, probably because of different atoms hav- 
ing different atomic weights. We were able to fit the data to the following formula which 
contains a linear combination of two Debye integrals. 
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CAT) = 9rNksJ2a {£) jj^^.dx (3) 

1-1,2 Q 

Here r is the number of atoms per formula unit, is a Debyc temperature and = 9^j^/T 
. The obtained parameters are Ci = 0.20 ± 0.01, 9}) = 183 ± 2 K, C2 = 0.41 ± 0.01, and 
= 456 ± 6 K. Further we also fit the Cp{T) of BiZnsPOe data to the low-T Debye 
expression {^NrkB (^^) ) in the T-range from 2 to 20 K. The yielded Od is (325 ± 2) 
K and is intermediate to 9]^ and 9%. In the same T— range, we go ahead and fit the Cp{T) 
data of BiCu2P06 using the following equation 



(4) 

Here the yielded value oi9D — 330 ± 2 from the first part is nearly the same as the value 
of 9d of BiZn2P06. Hence we used Cp{T) of BiZn2P06 as the lattice contribution for the 
BiCu2P06 (as mentioned before) and the magnetic part {x = 0.0) is plotted in Fig. 6. 

Exponential behavior, a broad maximum, and a lack of LRO are the indications of spin-gap 
systems. The second term is the magnetic specific heat for two-leg spin ladders. A fit of 
our data to the above Equation yields A = 43 ± 2 K (in good agreement with the 45 K value 
from low-T susceptibihty analysis) and 7/A = 6± 1. The magnetic entropy is estimated 
by integrating the data ^^21 ys T. The normahzed (S^) is plotted as a function of 
reduced temperature T*(= kBT/Ji) in the inset of Fig. 6. The value of S^^ at T* = 0.5 is 0.2 
which is 35% of the value of 0.56 for an 5" = | uniform chain which is non-frustrated.^^ The 
quenching of magnetic entropy at low temperature is a hallmark of frustrated spin systems. 



-A^ 



B. Zn-doped BiCu2P06 



The magnetic susceptibilities of Zn^+ {S = 0) doped at Cu^+ site in the BiCu2P06 
samples i.e., BiCu2(i-s)Zn2sP06, (0 < x < 0.5), were measured in an applied field of 5000 
Oe as shown in Fig. 5(a). The susceptibility at low temperatures increases with the doping 
concentration x. Introduction of Zn^+ at Cu^"*" in the ladder creates spin- vacancies, which 
are expected to induce localized magnetic moments on the neighboring Cu sites. The effect 
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of induced moments has certainly been observed by us in NMR and //SR measurements.^^ 
For the case of an isolated S — 1/2 HAF 2-leg ladder, the effective magnetic moment per 
vacancy is expected to be exactly that due to an 5" = 1/2 entity. For x > 0.02, it is clear 
that the low— T (T < 6 K) behavior is not Curie-like (there is a deviation from linearity in 
the log-log plot of Fig. 5(a)). We analyse the low-T behaviour in the susceptibility in the 
following manner. There are both intrinsic (natural defects/chain breaks, impurity- induced 
effects, spin-gap) and extrinsic (impurity phase) contributions to the susceptibility in all the 
samples. As argued, the spin susceptibility due only to the spin-gap (xnmr) is obtained 
from the ^^P NMR shift. Next, we subtract Xnmr from the measured susceptibility x- 
This yields Ximp(a^, T) which is the combination of a T-independent part and a T-dependent 
part due to defects, dopants, and extrinsic impurities and is shown in Fig. 5(b). As seen, 
the undoped sample has a weakly T-dependent contribution. We then plot xfmp(^ +^) 
divided by the Zn concentration, as a function of T in the inset of Fig. 5(b). The impurity 
susceptibility is xfmp = Ximpix) ~ Ximp — 0) ~ Xo where Xo is a T-independent term. Here, 
9 is taken equal to Tg which is the temperature below which a frozen magnetic state was 
evidenced from /xSR measurements. For a pure Curie- Weiss behaviour, Xknp(^ + ^) should 
be T-indcpcndent giving the value of the Curie constant. Our result is seen to be not quite 
T-independent (the nature of the data are not much dependent for 9 between and Tg). 
However, Ximp increases linearly with Zn content. Also, it is seen that the impurity-induced 
susceptibility is smaller (consistent with an effective paramagnetic moment of 5" ~ 0.3) than 
the case where each Zn induces an effective S — 1/2 moment (the dashed hue in the inset 
of Fig. 5(b)). This might be related to the presence of frustrating interactions in the 
compound. While the details of ^^P NMR measurements in pure and doped BiCu2P06 are 
discussed in a separate paper, suffice to say at this point that the data are consistent with 
the nominal Zn/Ni content in the ladder. The position of the broad maximum in x{T) does 
not appear to change with Zn content as shown in Fig. 5(a). Surprisingly even for x — 0.5, 
a broad maximum is seen. No difference in the zero-field-cooled (ZFC) and field-cooled 
(FC) susceptibilities was seen at low— T. However, the deviation from the Curie-behavior 
seems to suggest transition to a new magnetic state. These results are different from those 
of SrCu203, in which cusp like anomalies were found in x{T)^^ 

Quite generally, the replacement of Cu^"*" ion with Zn^"*" in a spin-liquid with a sin- 
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glet ground state like in HAF two- leg spin ladders induces an effective spin-| moment. 
However, a deviation from the free, spin-| moment was predicted in 5" = doped, one- 
dimensional, dimerized, frustrated (Majumdar-Ghosh) chains.^^ Considering that BiCu2P06 
can be thought of as consisting, essentially, of coupled Majumdar-Ghosh chains, predictions 
by Normand and Mila might be valid here. The reduced Curie constants (compared to the 
S — 1/2 value) observed by us in Zn-doped samples might be an evidence of this. However, 
there may be other explanations. For higher doping concentrations, interaction between 
the induced moments result in deviations in the x(T') from a Curie law. The origin of the 
disordered magnetic state which has been observed may be connected to random exchange 
interactions between impurity induced local magnetic moments on Cu sites on ladders in 
different 6c-planes. 

In order to better understand the properties of Zn doped BiCu2P06, we have also per- 
formed the specific heat Cp{T) measurements for x < 0.30. After subtracting the lattice 
contribution from the data for the non- magnetic analog BiZn2P06, the obtained magnetic 
contribution C^/T versus T is plotted in Fig. 6. Here we found that even for small doping 
level X = 0.01, an upturn is seen at low-T. This is possibly indicative of an onset of magnetic 
freezing. For higher x, broad anomaly is seen around 4 K. These suggest the loss of the 
spin-gap and a freezing of impurity-induced moments. Due to the small magnetic heat 
capacity and the weak anomaly, we don't have good accuracy in determining the freezing 
temperature. However, the anomalies in C^/T vs T are in the same range as found by our 
/iSR measurements on the same samples. 

C. Ni-doped BiCu2P06 

The x(T) of Ni^+ {S = 1) doped BiCu2P06 were measured in the T— range 2 to 300 K 
and are shown in Fig. 7(a). At low— T, an enhancement of susceptibility with Ni content {y) 
was observed and for higher contents (y > 0.15), even the high— T susceptibility increased 
with y. At low doping levels, we again carried out the analysis as for Zn doping. The plot 
of XimpC^ + ^) Ni-doped samples is shown in Fig. 7(b). A Ni dopant is expected to 
give rise to two contributions; (i) a staggered magnetisation on the neighbouring Cu sites 
as for the Zn and (ii) the magnetic moment of the Ni itself. While the "effective" Curie 
constant induced per percent of Ni (about 0.006 cm^K/mole) is clearly greater than that 



10 



for Zn, the value is smaller than that expected for 1% of paramagnetic S — 1 entities. 
This indicates that the Ni^+ impurity moments may be coupled antiferromagnetically with 
induced moments on neighboring Cu sites. However, it is also possible that there is a 
ferromagnetic coupling between the Ni and the neighbouring Cu. An agreement with our 
experimental results can be achieved if there is as well a frustration related screening which 
decreases the size of the paramagnetic cloud comprising of the Ni moment and the induced 
moments on the neighbouring Cu. Similar to Zn doped samples, at low temperatures x(T') 
deviates significantly from a Curie behavior. In addition, there are differences observed in 
ZFC and FC x{T) data measured in an apphed field 100 Oe as shown in inset of Fig. 7. 
Specific heat data arc further measured to understand these anomalies. Sharp cusp like 
anomalies are found with Ni doping as shown in Fig. 8. The temperatures at which these 
anomalies are present are similar to those in x(T) data. The exponential T— dependence and 
the broad maximum present in the Cm{T) of undoped composition were lost on Ni-doping 
indicating a loss of the spin-gap. The cusps in the heat capacity are likely associated 
with spin- freezing (similar features have been noticed in Zn- and Ni-doped SrCu203 [Ref.^^] 
and CuGeOs [Ref.^^]). Note that the susceptibility or heat capacity anomalies are not 
sharply defined in our data (in fact, if one looks at Cm vs T, only a broad feature is seen). 
Keeping this in mind, in our /xSR measurements (where the freezing temperature is obtained 
precisely) on the same samples, we observed spin-freezing at similar temperatures as the 
heat capacity anomalies. 

D. Ca/Pb-doped BiCuaPOe 

The replacement of Ca^+/Pb^+ in place of Bi^+ should generate holes. The x of hole- 
doped samples Bii_2(Pb/Ca)2Cu2P06 (0 < -2 < 0.15) are shown in Fig. 9(a). The x 
increases with doping concentration {z) and no anomalies were found at low-T, in contrast 
with Zn/Ni doped samples. Specific heat data are also consistent with this, as shown in 
the inset of Fig. 9. There is no change in the position of broad maximum with doping. To 
obtain Curie constant of Ca doped samples, we again plotted Xi^pT as a function of T in 
Fig. 9b. For clarity, we have not shown the data of Pb doped samples in Fig. 9(b). The 
effective Curie constant per percent of doped Ca/Pb amounts to less than one-fourth of that 
due to an S — ^ moment for each hole, which is even smaller than the value obtained for 
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Zn doped samples. This might be expected since Pb or Ca substitution at the Bi site does 
not alter the ladder geometry but is rather expected to dope holes. It is not known what 
type of spin-state these holes (which are evidently localised) form due to hybridisation of the 

Pb/Ca orbitals with the neighbouring orbitals. The Curie term might then be partly due 
to oxygen non-stoichiometry and partly due to the localised holes, apart from any extrinsic 
effects. 



V. CONCLUSION 



In summary, we have given a detailed account of preparation and properties of Ni-, 
Zn-, Pb-, and Ca-doped BiCu2P06. While Zn has complete solid solubility, Ni can be 
substituted upto a concentration of about 20 % and Ca up to about 15%. The effective 
magnetic coordination number for BiCu2P06 was found to be about 7 from an analysis 
of the susceptibility data. The heat capacity for BiZn2P06 is well fit in the T-range 2 
K-lOO K by a superposition of two Debye functions with Debye temperatures of about 
180 K and 450 K. The magnetic entropy of BiCu2P06 is quenched at low temperatures, 
presumably due to frustration. Both, Ni- and Zn-doping in BiCu2P06 are found to give rise 
to static, disordered magnetism at low temperatures. While clear anomalies were found in 
the temperature dependence of magnetic susceptibility and magnetic heat capacity of the 
Ni-doped samples, the Zn-doped samples exhibited merely a deviation from Curie behaviour 
at low temperature while the heat capacity showed a weak anomaly. The impurity-induced 
moment is found to be reduced from the S = 1/2 value in Zn-doped BiCu2P06 which might 
be due to frustration effects. On the other hand, Ca- and Pb-doped samples gave rise to 
only weak Curie terms and a nearly unchanged magnetic heat capacity C„i{T) which implies 
that the magnetic ladders/planes are not disturbed by the doping at Bi site and any carriers 
generated localise at low temperatures. 
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Figure Captions 

FIG. 1 (Color Online) (a)The crystallographic 6c— plane of BiCu2P06 is shown and it can 
be seen that coupled two-leg zigzag ladders formed by Cu (blue) & O (red) atoms and Bi 
(green) atoms are positioned in between. The various major exchange interactions between 
Cu atoms are shown in the schematic Fig. 1(b). 
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FIG. 2 (Color Online) The x-ray diffraction pattern of BiCu2P06 is shown along with the 
Rietveldt refinement. The black open circles indicate the experimental xrd data, the solid red 
line is calculated based on the refinement, green vertical short lines are the Bragg positions 
and blue line is the difference of measured and calculated data. The refinement parameters 
are Rp= 14.0, R^p = 11.0 and = 4.2. (b) The xrd patterns of BiCu2(i_a:)Zna,P06 {x = 0.0, 
0.5 & 1.0) with hkl-indices are shown, (c) The lattice constants a, b, and c as a function of 
Zn doping content (x) are shown in Fig. 2c. 

FIG. 3 (Color Online) (a)Magnetic susceptibility x{T) ^{M/H) vs. temperature T 
(open squares) for BiCu2P06 in an apphed field of 5000 Oe. Also plotted is XnmrIT") (open 
triangles). The solid line is a fit to Equation (2) of XiNMR(^)- As shown in the inset Fig. 
3(a), dashed line is a fit to the Equation (1) as explained in the text. In Fig. 3(b), XspinT 
is plotted as a function of T and inset shows magnetization (M) as a function of magnetic 
field (H) upto 90 kOe and a comparison with the Brillouin function at 350 K. The arrow 
mark indicates the value corresponding to the Curie constant (C) for free spin-|. 

FIG. 4 (Color Onhne) (a) The specific heat Cp as a function of T for BiCu2P06 (blue 
open up-triangles) and BiZn2P06 (black open down-triangles) is shown and fitted with the 
equations mentioned in the text. Inset shows the data of Cp as a function of for BiZn2P06. 
(b) The magnetic contribution to the specific heat of BiCu2P06 (left axis) is plotted vs. the 
normalised temperature k^T / Ji. Also shown is the magnetic entropy Sm (right axis) vs. 

FIG. 5(a) (Color Onhne) Magnetic susceptibility vs. temperature x(T') for 
BiCu2(i-x)Zn22;P06, (0 < x < 0.5) in the range 2 K to 300 K measured in the applied field 
5000 Oe. (b) A plot of the non-spin-gap part of the susceptibility Ximpfor BiCu2(i-2:)Zn22;P06 
(i.e., only from impurity-induced effects and from extrinsic impurities) is shown as a function 
of T. The inset shows Ximp(^ + ^) P^r percent Zn as a function of T. The dashed horizontal 
line represents the value corresponding to an 5" = 1/2 induced moment per percent Zn. 

FIG. 6 (Color Onhne) Magnetic specific heat divided by temperature Cm{T)/T vs. T for 
BiCu2(i-x)Zn2^P06, (0 < a; < 0.3) in the range 2 K to 50 K. 

FIG. 7(a) (Color Onhne) The x{T) of BiCu2(i_j/)Ni2j^P06 (0 < y < 0.2), measured in a 
magnetic field of 5000 Oe and plotted on a log-log scale. The ZFC and FC susceptibility 
x{T) measured in 100 Oe is shown in the inset. The arrow marks at bifurcation points 
are indicative of the freezing temperature, (b) A plot of the non-spin-gap part of the 
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susceptibility Ximpfor BiCu2(i-a;)Ni2a;P06 (i.e., only from impurity-induced effects and from 
extrinsic impurities) is shown as a function of T. The inset shows X^^{T -\-9) per percent 
Ni as a function of T. The dashed horizontal hne represents the value corresponding to an 
S = 1 induced moment per percent Ni. 

FIG. 8 (Color Online) Magnetic specific heat divided by temperature Cm{T)/T vs. T for 
BiCu2(i-2/)Ni2j/P06, {0 < y < 0.2), in the range 1.8 to 50 K (plotted on a semi-log scale). 
Anomalies are seen (arrow marks) which might be associated with magnetic freezing. 

FIG. 9 (Color Online) (a) The x{T) of Bii_^(Ca/Pb);,Cu2P06 {0 < z < 0.15) is shown. 
Inset shows CmiT) for pure and Pb-doped (z = 0.07) samples, (b) A plot of the non-spin-gap 
part of the susceptibility Ximpfor Bii_2(Ca/Pb)^Cu2P06 (i.e., only from impurity-induced 
effects and from extrinsic impurities) is shown as a function of T. The inset shows Xi^^p^ 
per percent Ca as a function of T. 
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